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ABSTRACT 


SMITH, WAYNE FARRIOR. Investigation of the Statistics of Pure Tone 
Sound Power Injection from Low Frequency, Finite Sized Sources in a 
Reverberant Room, (Under the direction of JAMES RONALD BAILEY and 
FRANKLIN DELANO HART). 

The measurement of sound power injection into a reverberant room 
is subject to statistical uncertainties owing to large fluctuations of 
sound pressure level with variations in source frequency and source 
and/or receiver position. Especially with regard to low frequency pure 
tone sources, point source theory indicates a requirement of many micro- 
phone and source positions if radiated power is to be accurately measured. 

An investigation has been conducted to determine the effect of 
finite source size on the power statistics in a reverberant room for pure 
tone excitation. Theoretical results indicate that the standard devia- 
tion of low frequency, pure tone finite sources is always less than that 
predicted by point source theory and considerably less when the source 
dimension approaches one-half an acoustic wavelength or greater. 

A supporting experimental study was conducted utilizing an eight 
inch loudspeaker and a 30 inch loudspeaker at eleven source positions. 

The resulting standard deviation of sound power output of the smaller 
speaker is in excellent agreement with both the derived finite source 
theory and existing point source theory, if the theoretical data is 
adjusted to account for experimental incomplete spatial averaging. 

However, the standard deviation of sound power output of the larger 
speaker is measurably lower than point source theory indicates, but is 
in good agreement with the finite source theory developed in this work. 
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1. INTRODUCTION 


Noise has long been recognized as a cause of annoyance and physio- 
logical damage to the human body. Before 1960, sporatic attempts were 
made to study and reduce this noise pollution problem, but no overall 
trend in this direction was manifest. However, recent public pressure 
and interest in noise has prompted considerable noise reduction re- 
search and legislative action. This new found interest in conquering 
the noise problem is exemplified by a special awareness toward the 
elimination of occupational deafness due to intense noise exposure. 

The Williams-Steiger Occupational Safety and Health Act of 1970 (25) 
has established a list of noise criteria which must be met by 
industries engaged in interstate commerce. The industries involved 
must comply with the established requirements or be subject to severe 
economic repercussions. 

One important compliance criterion of the Occupational Safety 
and Health Act (OSHA) involves the initiation of an engineering program 
to study and to reduce sound levels at their sources. This engineering 
study would, of course, include either retrofitting or redesign of 
existing noisy machinery or the purchase of new machinery designed to 
meet compliance standards. In either case, the ultimate objective 
would be to reduce the sound pressure level within the industrial 
environment at locations occupied by working personnel. Several indus- 
tries have been designated as target areas in the enforcement of the 
OSHA. These include longshoremen, roofing and sheet metal, meat 
products, mobile home manufacturing, and lumber and wood products. 
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The potential economic impact upon these industries of the strict 
enforcement of this act cannot be overemphasized. 

But the problems associated with any program designed to reduce 
noise levels involves the ability to evaluate accurately the acoustic 
benefit of retrofit and redesign or the ability to purchase new equip- 
ment with confidence that noise levels in the acoustic environment 
will truly be improved. Modifications made to existing machinery 
mounted in its normal operating position in the acoustic environment 
presents little, if any, problem, Comparative sound level measurements 
taken at normal personnel locations will indicate any beneficial 
changes. However, for any redesign done remote from the actual machine 
location or with regard to prospective machine purchases, sound pres- 
sure level alone is not a satisfactory measure of noise source strength. 
This is due to the fact that sound pressure level varies with the 
distance between the source and receiver and with the environment for 
both. The acoustic power and directivity of a source are the funda- 
mental quantities needed to describe acoustic performance. 

The acoustic power of a sound source is a measure of the total 
energy per unit time radiated in all directions, stated as a function 
of frequency. The acoustic directivity of a source is a measure of the 
change in power radiation magnitude with changes in angular position, 
also stated as a function of source frequency. Given the acoustic 
power and directivity, it is possible to make an engineering estimate 
of the resultant sound pressure levels generated by a source in almost 
any acoustic environment. Moreover, the acoustic power is often used 
by itself to compare the acoustic merits of different source designs. 
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In addition to the obvious beneficial application to machinery 
noise reduction research, acoustic power measurements could also play 
a big part in the establishment of product noise codes in the near 
future. It is anticipated that, due to annoyance considerations, such 
products as household appliances, lawnmowers, air conditioners, and 
motor vehicles will be required to meet noise emission standards. 

Sound power measurements could play a significant role in the estab- 
lishment of such noise standards. Thus, the need to understand and to 
develop an accurate and economical method of determining the sound 
power output of any given source is quite evident. 

In order to measure sound power with a high degree of confidence, 
special acoustical environments are required. One type of environment 
which is acceptable is termed anechoic. An anechoic facility is 
constructed by treating all the bounding surfaces of a room with sound 
absorbing wedges so that approximately 99 percent of the incident 
acoustic energy is absorbed within a typical frequency range from 
150 Hz to 10,000 Hz. The source under test is usually placed on a 
suspended floor within the facility and measurements of sound pressure 
level are made at specific locations on an imaginary hemisphere sur- 
rounding the source. Due to the free field environment existing in 
such a room, the sound level measurements enable one to accurately 
determine the free field sound power output of the source, in addition 
to its directivity characteristics. Such directivity is extremely 
important in the prediction of the noise level at any given location 
in an actual industrial environment. Obviously, the complete and 
accurate description of the sound power characteristics of a compli- 
cated source in an anechoic facility generally requires the 
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assimilation of a considerable quantity of data. Thus, in terms of 
man-hours, the test effort involved is relatively high. Additionally, 
the design and construction of an anechoic chamber requires a rela- 
tively high monetary investment. As a result, only a few anechoic 
facilities are in existence in the United States and these are by no 
means readily available for the volume of sound power measurements 
anticipated in the near future. 

The converse of the anechoic environment is the reverberant envi- 
ronment. A reverberation chamber is designed so that the bounding 
surfaces absorb very little acoustic energy, resulting in a highly 
reflective, diffuse sound field. Generally speaking, sound power 
measurements in a diffuse sound field require less test effort due to 
the fewer number of microphone readings and the less accurate microphone 
placement required. For sources radiating essentially broad band noise, 
the diffuse field sound reflections average out in the room and the 
average of only a few microphone readings is well representative of 
the spatial average pressure in the room. Sources which exhibit 
strong pure tones or narrow bands of noise, however, can generate 
standing waves, thereby resulting in large fluctuations in sound pres- 
sure with variations in source frequency as well as source and/or 
receiver positions. Associated with these fluctuations and variations 
is a certain amount of statistical uncertainty which must be fully 
understood in detail before pure tone measurements can be undertaken 
with confidence. Thus, the reverberation room technique, which was 
first thought to be an effective tool for acoustic research, has 



5 

developed into a research obstacle which must first be understood 
before it can be effectively used for pure tone measurement. 

The state-of-the-art for reverberant rooms is fairly well estab- 
lished for broad-band, narrow-band, and pure tone excitations if the 
source is assumed to exhibit monopole type radiation characteristics 
and to be located at a point in the reverberant field. However, the 
majority of real sources are not located at a point but are distributed 
over some finite area. Therefore, this study is designed to investigate 
the effects of finite source size on the low frequency variance of 
acoustic power for a monopole type pure tone source in a reverberant 
room. The theory developed would be equally applicable to high fre- 
quency variance; but, as pointed out by Baade\ the total variance at 
high frequency is relatively small so that any change due to finite 
source size would be of little practical importance. 

The presentation of this study begins in Chapter 2 with a review 
of the present state of the art in reverberant room power measurements. 
In this literature review, special emphasis is placed on reverberant 
rooms excited by pure tone sources. 

In Chapter 3, the output power statistics for a pure tone point 
source are reviewed. This theory is then extended in Chapter 4 to 
include finite source dimensions. The analysis of the finite source 
model is discussed in detail with the results expressed in terms of 
corresponding point source theory. 


This remark was made to the Acoustical Society of America as part 
of the introduction and discussion for the Symposium on Sound Power 
Measurements in Reverberant Rooms, November 22, 1968. 
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Chapter 5 provides a complete description of the reverberation 
facility which was utilized in obtaining experimental results in sup- 
port of the theoretical derivations. This description includes both 
theoretical and experimental calculations of room acoustic properties, 
as well as a thorough description of the experimental setup. Special 
emphasis is placed on a description of the microphone traversing 
apparatus which was specially designed and constructed for use in this 
study. Finally, the statistical analysis of the data utilizing a two 
way analysis of variance is discussed in detail and the corresponding 
results are presented. 

In Chapter 6, the experimental results are compared to both point 
source theory and the derived finite source theory. Special emphasis 
is placed upon the necessary adjustment of the theoretical results to 
account for the experimentally incomplete spatial averaging. 

Finally, in Chapter 7, the significant findings and contributions 
of this study are summarized and conclusions are drawn. These con- 
clusions include suggestions for expanding this study, as well as a 
brief discussion of some other problems, yet unsolved, associated with 
reverberant rooms. 
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2. REVIEW OF THE LITERATURE 

The study of reverberation environments and the properties of 
such environments which made them suitable for acoustic power measure- 
ments is only a small part of a larger area of research known as room 
acoustics. Although the acoustics of enclosed spaces has been of 
obvious importance since the first auditoriums and theaters were built, 
it was late in the nineteenth century before the maze of superstition 
surrounding the subject began to be replaced by scientific explanation. 
But major breakthroughs were still slow in coming; in fact, it was not 
until the last two decades that a basic understanding of reverberant 
field power measurements began to be developed. 

The first major development concerning the acoustic properties 
of large enclosures was made by Sabine (19) in 1895. Through extensive 
experimental studies, Sabine arrived at a mathematical relationship 
between the reverberation characteristics of an enclosure, its volume, 
and the amount of absorption present. This relationship is given by 


0.049V 

Sa 


( 2 . 1 ) 


where T is the reverberation time of the enclosure in seconds, V is the 
room volume in cubic feet, S is the room surface area in square feet, 
and a is the average absorption coefficient of the room. His defini- 
tion of reverberation time as the time in seconds required for a 60 dB 
decrease in sound intensity is still the single, most important param- 
eter in judging the quality of an acoustic environment. Theoretical 
derivations of (2.1) are usually based on geometric acoustics in which 
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sound is assumed to travel in rays (17). At each encounter with a 
boundary, the rays are partially absorbed and partially reflected. 
Although this theory often oversimplifies the actual behavior of waves 
in rooms, it leads to valid results under certain specific conditions. 
If the absorptive capability of the room is high, the number of wall 
reflections is small, the reverberation time is small, and the 
ultimate sound intensity at any point is quickly reached. Such rooms 
deviate greatly from Sabine's theory. On the other hand, if the 
absorption is low, the intensity growth rate is small and the rever- 
beration time is large. Rooms of this type are known as live or 
reverberation rooms, and application of Sabine's theory produces 
accurate results. 

For sound sources started in a live room such as described above, 
many reflections at the walls produce a sound energy distribution that 
becomes more and more uniform with time (16). Ultimately, this energy 
density may be assumed to be uniform, with completely random directions 
of flow. Such a condition, known as diffuse, is extremely important 
in room acoustics due to its relative ease of description using normal 
mode and statistical techniques. 

Large rooms under diffuse conditions can be used to accurately 
measure the absorption and transmission loss of acoustical materials 
and the sound power output of acoustical sources. The problem, however, 
which arose early in reverberation room research, and a problem which 
still exists to some extent today, is the establishment of a valid 
criterion to determine whether the sound field in an enclosure can be 
assumed to be diffuse (21). A large reverberation time would tend to 
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indicate diffuse conditions, but the intimate relationship among room 
volume, excitation frequency, and the amount and location of absorption 
prevents the establishment of any limiting value for an acceptable rever- 
beration time. 

The simplest method to date for determining the existence of diffuse 
conditions involves the calculation of Schroeder’s "large room frequency" 
(20) , given by 


4 1 1/2 

f s - 1.19 x 10 <£) 


( 2 . 2 ) 


where T and V are as defined for (2,1). This quantity describes the low 
frequency cutoff above which the sound in a particular room may be approx- 
imated as diffuse. Experimental results, however, have shown this number 
to be quite conservative, in that data taken one octave or more lower 
has exhibited satisfactory diffuse characteristics (7) . 

To help alleviate the obvious problems which exist concerning the 
diffuse field approximation, ASHRAE Standard 36-62 is often utilized 
(1,22). But this qualification procedure lacks consideration of recent 
developments and is, consequently, rather incomplete as far as pure tone 
testing is concerned. An improved test procedure, therefore, is presently 
under draft which hopefully will define completely both the usable fre- 
quency range of any reverberation facility and the test procedure for 
pure tone and narrow-band testing. 

With the importance of diffusiveness firmly established, several 
other parameters and their effect on reverberant room power measurements 
need to be discussed. A number of studies (6, 7, 12, 14, 20, 26) have 
been carried out on the fluctuations in room sound pressure when the 



frequency of the source is varied or when the position of the source 
or receiver is varied. These fluctuations give rise to statistical 
uncertainties which must be resolved and understood before accurate 
measurements can be performed. For broad-band noise, theory and experi- 
ment indicate very little uncertainty, except in cases of very low room 
absorption (2). This is because the spatial fluctuation of pressure 
throughout the room is small and the change in the reaction of the room 
environment on the source is quite small. However, problems arise in 
trying to measure the sound power of sources exhibiting discrete fre- 
quencies. Much of the current research in room acoustics has been aimed 
at understanding and solving the problems of pure tone measurements. 

The statistics of source irregularity (frequency irregularity) , 
originally studied by Bolt and Roop (6), was revised and extended by 
Schroeder (20) . He found that if the large room frequency was satis- 
fied, the standard deviation from the mean of room sound pressure level 
due to changes in source frequency was 5.57 dB. Furthermore, it was 
found that the lesser the degree of room symmetry, the greater was the 
degree of sound diffusiveness. The room irregularities decrease the 
concentrated areas of degenerate modes in the room frequency response. 

In 1959, Doak (7) investigated the fluctuations in sound pressure 
level when the receiver position is varied - space irregularity. Using 
both theoretical and experimental techniques, the standard deviation of 
sound pressure level from its mean for a pure tone source was found to 
be the same for variations in receiver position as it is for variations 
in source frequency, namely 5.57 dB. This value, however, is only valid 
for receiver positions far from the source and beyond one-quarter of a 
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wavelength from the room walls and for frequencies near or above the 
large room frequency. 

In several later publications, Lubman expanded upon Doak's work to 
better describe the fluctuations in pure tone pressure due to receiver 
position. In 1968 (12) , he showed that the standard deviation of pure 
tone pressure normalized to the mean has a value of unity. This was 
followed in 1969 (13) by a discussion of the reduction in normalized 
standard deviation by a ratio of vf if N independent samples of the 
sound field are averaged. These independent samples can either be taken 
discretely at half-wavelength intervals or be obtained by averaging over 
a straight line traverse. The equivalent number of independent samples 
in a traverse is given by N = 1 + 2X/X, where X is the length of the 
traverse and X is the acoustic wavelength. Experimental investigations 
have shown that traverse averaging is as good or better than discrete 

2 

point averaging, in addition to frequently being a time saving procedure. 

The final statistical variable which must be considered in reverber- 
ant field power measurements is that due to variations in source position. 
The acoustic power output of any given source is a function of the acous- 
tic impedance presented to it by the surrounding medium. The modal 
patterns set up in a reverberation room may increase or decrease this 
impedance depending upon source location and upon source frequency in 
relation to room modal frequencies. Lyon (14) determined the theoretical 
variance in average pure tone sound power with changes in source position 

2 

This idea is from the paper "Experimental Evaluation of Microphone 
Traverses in a Reverberant Room" presented to the Acoustical Society of 
America by C. E. Ebbing and D. J. Ingalls on November 7, 1969. 
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for a monopole type point source. His procedure and results will be 
discussed briefly in the next chapter. Computer studies (9, 15) have 
confirmed Lyon's results and have also shown that the greatest contribu- 
tion to the variance comes from source positions within one-quarter of 

a wavelength of the chamber walls, floor, and ceiling. 

3 

Ebbing (9) and Tichy have obtained lower variances by installing 
a rotating diffuser in the reverberation room to continuously shift the 
modal patterns. By use of a properly designed diffuser, two and even 
threefold decreases in variance have been obtained. 

With the major sources of error in reverberant room power measure- 
ments defined, the total variance is simply the sum of all individual 
variances (2). Thus, to predict the total uncertainty in a given power 
measurement, the variances due to all significant sources of error must 
be thoroughly understood. Once understood, of course, this total uncer- 
tainty can be decreased to an acceptable level by the choice of suffi- 
cient measurement locations . 


3 

J. Tichy presented the paper "The Effect of Rotating Vanes on the 
Sound Field in Reverberation Chambers at Single Frequencies" to the 
Acoustical Society of America on November 3, 1970. 
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3. POINT SOURCE THEORY 


The statistics of power injection associated with a monopole type 
point source in a diffuse field already exist in the literature. This 
theory is briefly reviewed in this chapter due to the fact that the 
finite source development in this study is derived as an extension of 
point source theory. 

From Maling's (15) development of a harmonic point source in a 

reverberation room, the sound pressure at any point (x,y,z) due to a 

monopole type source located at point (x , y , z ) is given by 

s s s 


P(x,y,z,t) 


8ipck(Qoe j ~ Uit ) r ^m(x,y,z)^m( X s, y s ,Z s) 

V 1 2 2 

V m (lc -k ) + 2igkS 

m — — — 


(3.1) 


where 

= modal amplitude function for a point source, 
m 

Qoe = source volume velocity, 

oi = source angular frequency, 

t = time, 

k,k m = excitation and modal wave numbers, respectively , 

g = damping coefficient, 

p = density of air, 

c = speed of sound in air, 

V = volume of room, 

S = surface area of room. 
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Evaluation of equation (3.1) at the surface of the source and 
division by the source volume velocity yields the input acoustic imped- 
ance of the source. This acoustic impedance indicates the reaction of 
the acoustic medium upon the vibratory motion of the source. The real 
part of this impedance, the acoustic resistance R, is associated with 
the dissipation of acoustic energy. Due to the fact that the acoustic 
power radiated by a source is equal to the rate at which work is done 
against the acoustic resistance, this resistance governs the sound power 
output of the source. Therefore, determination of the variance of in- 
put acoustic resistance due to changes in source position is analogous 
to determining the actual variance in source power. The expression 
for the input acoustic resistance of a harmonic point source is 


R(x ,y ,z 
s’-'s* s 


) - 16pck 6S £ 


m 


, 2 , , 

♦m (WV 

(k 2 -k 2 ) 2 + ,2k$S 2 

111 ' T7 ' 


(3.2) 


Using equations adapted from Rice (18) , Lyon utilized an expression 
similar to equation (3.2) to obtain the mean acoustic resistance and 
the variance of acoustic resistance normalized to the mean for a pure 
tone point source in a diffuse field. The results are 

k 2 pc 

m R = — » (3.3) 


2 

g R = 2 2_iZ 1 

2 °N 16 M » 
“R 


(3.4) 
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where 

m = mean acoustic resistance, 

K 

2 

a = variance of acoustic resistance or 
K 

acoustic power, 

2 

a XT = normalized variance of acoustic 
N 

resistance or acoustic power, 

M = modal overlap of room. 


The modal overlap function M is given by 

M _ Trkn(k)n 
M 2 


(3.5) 


where n(k) is the wave number modal density of the room and n is the 
damping loss factor. 

It should be noted that the mean acoustic resistance in equation 
(3.3) is equal to the acoustic resistance of a point source in a free 
field. This means that the average acoustic power output of a point 
source in a diffuse field is equal to its power output in a free field. 

The normalized variance in equation (3.4) can be put in a more 
practical form by expressing the modal overlap in equation (3.5) in 
terms of measurable acoustic properties. The damping loss factor is 
assumed to be equal to (23) 



where f is the source frequency in Hertz. The modal density of a 
rectangular room in terms of the wave number is (3) 
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n(k) 



(3.7) 


Substituting equations (3.6) and (3.7) into equation (3.5) the 
following expression for M is obtained: 


M = 


4.47T 2 f 2 V 

c 3 T 


(3.8) 


This expression for M is then substituted into equation (3.4) to 
give an expression for the normalized variance in terms of measurable 
properties. The result is 


7.39 x 10 3 ,T n 1/2 
ct n = f "v 


(3.9) 


In order to get a better grasp of equation (3.9), assume the 

~ 3 

existence of a reverberation room of scale length L, so that V = L . 
Then, from equation (2.1) 

V 1/3 

T - 0.0081 f . (3.10) 

a 

Substitution of equation (3.10) into equation (3.9) yields 


665 1 

°N _l/2 1/3 

a fV 


(3.11) 


which is graphed in Figure 3.1. This graph clearly shows the effects 
of absorption coefficient, frequency, and room volume upon the standard 
deviation of pure tone sound power for two typical values of a. An 
increase in either average absorption coefficient, frequency, or room 
volume decreases the normalized standard deviation of power. 
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cr 

^ 0.9 

2 0.8 




f V l/3 x I0 3 Hz • ft. 


3.1 Uncertainty of sound power output at a single arbitrary source 
position for a pure tone point source 
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4. FINITE SOURCE THEORY 
4.1 Mean Acoustic Resistance 

The previous analysis for a point source can now be used in part to 
determine the statistics of power injection for a monopole type source 
of finite size. The finite source is modeled as a baffled square piston 
of dimension d oriented parallel to the y-z plane as shown in Figure 4.1. 
The location of this piston in the reverberation room is made with refer- 
ence to the piston's geometric center which is given the arbitrary desig- 
nation (x g , y , z g ) . Considering this finite source to be composed of 
a number of point sources in phase, the desired input acoustic resistance 
for the source model can be obtained by averaging the acoustic resistance 
of each constituent point source, as given by equation (3,2), over the 
finite source area. This averaging process is represented by 

R(x s’ y s’ Z s ) = T2 R<X s» y s» Z s ) dA S * (4a) 

d U 

s 


where R is the average acoustic resistance of the finite source model 

and A is the surface area, 
s 

In order to perform the integration indicated in equation (4.1), an 
expression is needed for the point source modal amplitude function ij> m in 
relation to the source location. For a rectangular room, this function 
is given by (24) 

ib = cos k x cos k y cos k z 
m x s y s z s 




(4.2) 
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4.1 Finite source model positioned in a reverberation room 
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where k , k , and k are the components of the modal wave number k in 
x y z m 

the x, y, and z directions, respectively. Substitution of equations (3,2) 
and (4.2) into equation (4.1), with the establishment of proper limits, 
yields 


a , 2._ cos k x^ 

«VV*.> ■ I i ^ J r 

m 


2 2 
d V 


(k 2 -k 2 ) 2 + (^MV 

m V 


y +d/2 z +d/2 
s s 


2 '2 ' ' 
cos k Y cos k Z d d 
y z y z 


(4.3) 


y s“ d/2 z s “ d/2 


where y and z are arbitrary variables of integration. The results of 
the integration are 




m d2 y2 (k 2 -k m 2 ) 2 + (~^-) 2 


Y 1 sin 2y 

2k 4k 

y y 


k (y +d/2) 


k (y -d/2) 

y s 


z 1 sin 2z 

2k 4k 

z z 


k (z +d/2) 
z s 


k (z -d/2) 
Z s 


or 


R( ws> 


- 1 


16pck gS 


cos k x 
x s 


" d2 v2 0c 2 -k m 2 ) 2 + (^i) : 
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y c +d/2 sin 2k (y +d/2) 

(~ + ) - (■ 


y -d/2 sin 2k (y -d/2) 

V— + 3L.J ) 


4k 


4k 


z +d/2 sin 2k ^z +d/2) 

(— + 5 — s 

^2 4k ' 


sin 2k z ^z g -d/2) 


z -d/2 

(A— + 


4k 


(4.4) 


In order to simplify equation (4.4) for easier manipulation, let 


. 16pck 0S 

2 * 
V 


(4.5) 


B - (k 2 -k 2 ) 2 + (2Mi)‘ 
m V 


(4.6) 


Reduction of equation (4.4) then results in 


R (x ,y ,z ) = l cos“ k x 


s"s' s' ,2 

m d B 


2 V x 
x s 


[ 

3 

[ 


d sin 2k y (y g +d/2) sin 2k^(y g -d/2) 


2 + 


4k 


4k 


. sin 2k (z +d/2) sin 2k (z -d/2) 

Q X “ ° Z S 

2 4k 4k 


(4.7) 


From a knowledge of trigonometric relationships. 


sin 2(p+q) - sin 2(p-q) = sin 2p cos 2q 


+ cos 2p sin 2q - sin 2p cos 2q 


+ cos 2p sin 2q 
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and 


cos 2p = 2 cos p - 1. 


Thus , 


2 

sin 2(p+q) - sin 2(p-q) = 4 cos p sin 2q - 2 sin 2q 


(4.8) 


Utilization of the relationship of equation (4.8) in equation (4.7) 
yields 


R(x ,y ,z ) = T —x— cos k x 
s’^s’ s L ,2_ x s 

m d B 


, 2 cos k y sin k d - sin k d 

! + ( ZlS. I 2— \ 

Z K 2k ' 


, 2 cos k z sin k d - sin k d 

d . , z s z z 

2 * < 


2k 


■ . 


(4.9) 


The bracketed quantity in equation (4,9) can finally be put into non- 
dimens ionalized form so that the average acoustic resistance of the 
finite source model can be expressed either in factor form by 


R(x ,y ,z ) = y ~ cos 2 k x 
s s s L B Vi 


m 


x s 


. k d - sin k d _ sin k d 

l ( . y_ — 2_) + ™ 2 1 v 2- 


k d 

y 


-2— ) + cos k y 


s k d 

y 


- k d - sin k d - sin k d 

i,z z x , 2 , z 

_ ( — ) + cos ^ 2 _ — 

2 kd zskd 


(4.10) 



23 


or in expanded form by 


E( W 2 s> ■ l f 


m 


. 0 k d - sin k d 

7 - COS k‘x (— * : -j ^— ) 

4 x s k d 


k d - sm k d .. „ 0 sin k d 

, z z s , 1 2 , . 2 . z 

( ; J ) + cos k X COS k Z - j 

kd 2 xs zskd 


k d - sin k d 


k d 

y 


) + 4 cos 2 k x cos 2 k y 
2 x s y-'s 


sin k d 

2L. 

k d 

y 


kd - sin k d 000 

( r — ; ) + cos k x cos k y cos k z 

kd xs ys zs 


sin k d sin k d 

L. L_ 

k d k d 

y z 


If we represent the bracketed quantity in equation (4.11) by the 


symbol <|> and recall the values of A and B from equations (4.5) and 
(4.6), the finite source acoustic resistance takes on the form 


(4.11) 


R(x .y ,z ) = 


s s 


16pck gS 

„2 


m 


m 


(k 2 -k 2 ) 2 + 
m 


, 2gkS , ' 
^ V ’ 


(4.12) 


which is very similar to the acoustic resistance for a point source given 

o 

by equation (3.2). The only difference is the replacement of ib for a 

m 

2 

point source by <j> must represent the squared modal amplitude function 
for the finite source model. 

It should be noted that in the limit as the finite source dimension 
d approaches zero, equation (4,12) reduces to 
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lim R(x ,y ,z ) = 
s J s s 

d+0 


16pck gS 


2 2 2 
cos k x cos k y cos k z 
x s y s z s 


m „ 2 , 2 S 2 , / 2gkS > 


(k -k T + (- 

m V 


16pck gS 
.,2 


m 


m (k 2 -k V + 

m V 


Thus, as the size of the finite source becomes very small, the average 
finite source acoustic resistance reduces to that for a point source. 

A graphical sketch of equation (4.12) would exhibit randomly spaced 

peaks as shown in Figure 4.2, resembling a sequence of pulses from a 

random signal generator. The peak amplitudes would depend upon the 

2 

modal amplitude function <f> m and upon the source location, which is 
assumed to be equally probable at all points in the room. Furthermore, 
it has been assumed that all modes in a given frequency interval have 
approximately the same damping coefficient 3. This is equivalent to 
saying that all the peaks or pulses in a given interval have nearly 
the same shape, and 3 is not modal dependent. 

The statistics of pulse sequences with Poisson separations have 
been studied at length by Rice (18) , thereby making expressions for the 
mean and variance of such sequences readily available. Therefore, for 
the entirety of this theoretical derivation, the modal peaks will be 
assumed to exhibit this type of statistical distribution. 

Using Rice's formula, as adapted for special use in acoustic en- 
vironments by Lyon, the mean value of the acoustic resistance in 



8k= WAVENUMBER SEPARATION 



4.2 Graph of acoustic resistance in a room 
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equation (4.12) is given by 


M = 16 P ck g S 
V 2 


J 


<d> > n (k )dk 

T m ' m m 

,.2 . 2.2 ,2k$S ‘ 

0 (k " k m } + ( "T~ } 


(4.13) 


2 

The quantity <<J» > represents the spatial average of the squared finite 

source modal amplitude function. This spatial average is performed over 
all possible source positions in the room volume and can be represented 
by 



1 

V 


L L L 
x y z 

t r r 





. 


oo o 



dx dy dz , 
s s s 


(4.14) 


where L , L , and L are the room dimensions in the three coordinate 
x y z 

directions. 

Recall from equation (4.11) that 



, o k d - sin k d . k d 

i k d y > <— 


sin k d 
z . 


k d 
z 


. 7 o sin k d k d - sin k d 

+ cos k x cos x z — : — (—^ — : — j 2L- ) 

2 xs zskd kd 


. 0 o sink d k d - sin k d 

.12, 2, y,z z. 

+ - 7 T cos k X cos k y — ' ( t — j ) 

2 xs y's k d kd ' 


2 2 2 sin k d sin k d 

+ cos k x cos k y cos k. z — ; — — . 

xs y s zs kd kd 


(4.15) 
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Again for simplicity in manipulation, make the following substitutions 
in equation (4.15): 


k d - sin k d 

D = Z_ 

k d 

y 


e = 


k d - sin k d 
z z 

k d 
z 


G = 


sin k d 

JL. 

k d 

y 


H = 


sin k d 
z_ 

k d 
z 


(4.16a) 


(4.16b) 


(4.16c) 


(4.16d) 


These substitutions then enable one to rewrite equation (4.14) as 


, 2 1 

< 4 > > = 77 

m V 


L L L 
x y z 

* r * 


o o o 


DE 2. 

7— cos k X 
4 x s 


, DH 2, 2. , EG 2. 2, 

+ tt— cos k x cos k z + -r— cos k x cos k y 
2 xs z s 2 xs ys 


2 2 

+ GH cos k x cos k z 
x s z s 


dx dy dz 
s s s 


The integration is now carried out to give 


«f> m 2> = -| [DE + DH + EG + HG] . (4.17) 

Referring to the expressions for D, E, G, and H from equations 
(4.16a), (4.16b), (4.16c), and (4.16d), respectively, it can be shown 
that in the limit as d approaches zero, 
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D and E 0 , (4.18a) 

G and H -* 1 , (4.18b) 


and therefore, 


lim <<p 
d-*0 


m 


2 


> 



0 


2 

As might be expected, this limiting value of <<j> m > is equal to the value 
2 

of <i|i > obtained by Lyon for a point source. Thus, the spatial average 
of the squared modal amplitude function for the finite source approaches 
that for a point source as the finite source dimension becomes very 
small. 

2 

This value of <<j> >, however, instead of being a function of the 

modal wave number k^, is a function of the components of the modal wave 
number in each coordinate direction [k^, k y , k z ], Substitution of equa- 
tion (4.17) into equation (4,13), therefore, does not result in an inte- 
grand amenable to direct integration with respect to the variable of 
integration k m » A means of transforming equation (4.17) to an expression 
depending solely upon k m and not its vector components is definitely 
needed. 

One means of performing this transformation is through the utiliza- 
tion of a wave number lattice diagram. The relationship 




suggests that each modal wave number may be represented as a vector in 

wave number space (16), having components k , k , and k . Each of these 

x y z 



29 


components can be expressed in terms of the properties of the reverber- 
ant enclosure by 


n it 
k 

x L 

x 


(4.19a) 


n ir 

k = JL. , 

y L y 


(4.19b) 


k = n ir/L , 
z z z 


(4.19c) 


where n , n , and n represent the actual mode numbers in the x, y, and 
x* y z 

z directions, respectively. A normal mode of vibration can, therefore, 

be represented as a point in wave number space, with the x component of 

this point being some integral multiple of tt/L^, the y component some 

multiple of tt/L^, and the z component some multiple of it/L z . Such a 

wave number lattice for a rectangular room is shown in Figure 4.3. All 

normal modes having a characteristic wave number equal to or less than 

k are included in an octant of this space of radius k . 
m m 

This point lattice representation of the normal modes of a room can 
prove quite useful for making two approximations necessary to this study. 
First of all, the lattice can be used to obtain an equation approxi- 
mating the number of modes existent in a particular frequency interval. 

Each lattice point is assumed to occupy a space of dimensions tt/L x , 

3 

ir/Ly, and tt/L z and of volume it /V. The number of modes N below a given 

wave number k can then be determined by dividing the volume of an 
m 

3 3 

octant of radius k [irk /6] by tt /V. This yields 


k 3 V 
m 


(4.20) 
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A still more useful quantity can be obtained by differentiating 

equation (4,20) with respect to k . This yields the number of normal 

modes contained in the wave number interval dk , centered at k . The 

m m 

result is 



k 2 V 
m 



(4.21) 


which should be recognized as the approximate modal density n(k m ) of a 

rectangular reverberation room [equation (3.7)]. An increase in this 

modal density, either by an increase in wave number [higher frequency] 

or by an increase in room volume, results in a higher degree of random 

distribution of sound energy. As discussed in Chapter 2, any elevation 

in the degree of random distribution serves to improve the diffuse 

approximation for the reverberant field. 

The second important use of the wave number lattice is as a means 

of expressing equation (4.17) as a function of k . This is accomplished 

m 

through recognition of the fact that the components of any point in the 
wave number lattice can be expressed in terms of k m by use of spherical 
coordinates. The representation of a typical point is shown in Figure 
4.4. The results as obtained from the figure are 

k = k sin 6 cos y » (4.22a) 

x in 

k y = k m sin 9 sin y , (4.22b) 

k = k cos 6 . 

z m 


(4.22c) 




Substitution of these coordinate relationships into equations (4.16a), 
(4.16b), and (4.16c) then results in 
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k d sin 0 sin y - sin(k d sin 0 sin y) 
m m 


k d sin 0 sin y 
m 


(4.23a) 


k d cos 0 - sin(k d cos 0) 
m m 

k d cos 0 
m 


(4.23b) 


sin (k d sin 0 sin y) 
m 

k d sin 0 sin y 
m 


(4.23c) 


sin (k d cos 0) 
m 

k d cos 0 
m 


(4.23d) 


After the insertion of equations (4.23a), (4.23b), (4.23c), and 

2 

(4.23d) into equation (4.17), the dependence of <<t> m > upon k^, k^, and 

k z is eliminated, but the dependence upon the angular position of k m in 

the wave number lattice has been introduced. This new dependence problem, 

2 

however, can be solved by averaging <<J> m > over all possible angular loca- 

4 2 

tions of k^ in the lattice. This is equivalent to stating that <<j> m > is 

averaged over the surface of an octant of radius k in the wave number 

m 

TT / 2 tt/2 



1 

k 2 tt/2 
m 




sin 0 d 0 d y 


o o 


(4.24) 


This averaging process was suggested by S. N. Yousri and F. J. 
Fahy, University of Southampton, England. 



34 


where k n/2 is the surface area of the octant and k sin 0 d 0 d y 
m m 

is the area of a typical differential element. 

2 

Because of the extremely complicated nature of <<j> m >, equation 

(4.24) does not lend itself to solution by direct integration methods. 

However, an approximate integration can easily be performed by treating 

the integral as a finite double sum and by using numerical solution 

techniques. A description of the digital computer numerical integration 

procedure which was devised and utilized is given in Appendix 9.1. 

The value of equation (4.24) is obtained as a function if k d for 

m 

values of k d from 0.25 to 25.00 at intervals of 0.25. The computer 
2 

results show <<j) > to be a constant equal to 1/8 for all values of k d. 

m m 

2 

It should be noted that this computer value of <«j> > is equal to the 

2 

value of <il> > for a point source, 
m 

Utilizing the result 

< 4 > 2 > = 1/8 
m 

in equation (4.13), the following integral for the mean value of acoustic 
resistance is obtained: 


2pck 6S 
"R - v 2 


n(k ) dk 
m m 


(k 2 -k 2 ) 2 + 

m V 


(4.25) 


The damping coefficient B in equation (4.25) is related to the 
wall damping loss factor n by 


B 


kqV 

2S 


(4.26) 
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If n is small [<0.1], which it must be for diffuse conditions to be 
approximated, then the integral in equation (4,25) peaks sharply about 
k^ = k, and the following reduction can be made: 


2pcgS 


n(k) dk 


( k 2 -k 2 > 2 + 

m V 


(4.27) 


Letting 


R = 


2 2 
k - k 

m 


and 


2 2 
k + k 

dR = -5™ dk 

, 2 m 


= 2 dk 

m 


then equation (4,27) transforms to 


m R 


PC 

2V 


n(k) 


r 


j- CO 


( 23 s /V) 

R 2 + C^f ) 2 


dR 


for small damping. Integration of equation (4.28) yields 


"R = 


pcn(k)fr 

2V 


(4.28) 


(4.29) 


Recalling the approximate value of the modal density of the room n(k) 
from equation (4.20), the final result for the mean acoustic resistance 
of the finite source model is 


pck 

“r = — 


(4.30) 
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This is exactly the same result as was obtained by Lyon for a point 
source [equation (3.3)]. Thus, the mean acoustic resistance of the 
finite source model in a reverberant room is equal to that for a point 
source in a reverberant room and also equal to that for a point source 
in the free field. ^ 

This result takes on a greater significance if recognition is 
made of the relationship between average acoustic power and mean acoustic 
resistance. This relationship is given by 
2 2 

W = S V HL , 

K. 

where 

W = average acoustic power, 

S = source radiating area, 

_2 

V = mean square source velocity. 

2 _ 2 

For monopole type radiation, S V is equivalent to the source volume 
velocity squared. Thus, for the case of equal volume velocities, a point 
source and the finite source model exhibit the same average power output 
in a diffuse field. Furthermore, this average power is equal to the 
power output of a point source or the finite source model in a free field. 
This is to be expected if the reverberant room technique is to be used 
for measurement of the classic free field acoustic power. 


5 2 

The dependence of the mean acoustic resistance upon k [equation 

(4.20)] for values of ka >> 1 has been questioned by P. K. Baade. How- 
ever, this result has no influence upon the other significant conclusions 
in this work. 
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4.2 Standard Deviation of Acoustic Resistance 
As mentioned in Section 4.1, the variance or standard deviation of 
pulse sequences with Poisson separations has been studied in detail by 
Rice. Thus, using an equation developed by Rice and adapted for acous- 
tical use by Lyon, the variance of the average acoustic resistance for 
a finite source [equation (4.12)] is given by 



, 16pck^gS > ^ 

(. 2 ' 

V 


<<(>_> n(k ) dk 
m mm 


2 2 


rn 2 . 2.2 ,2kgSv . 
o [ (k -k ) + ] 

m v 


(4.31) 


The quantity <<(> > represents the spatial average of the finite source 

modal amplitude function raised to the fourth power. Referring to 

equations (4.15), (4.16a), (4.16b), (4.16c), and (4.16d), the value of 

2 L 

<f> m can be squared and simplified to yield the following result for <j> : 


2 2 2 2 

, 4 D E 4 , , TTYL 4 , 4 , 

<p = ~ r cos k x 1 — cos k x cos k z 

m 16 xs4 xs zs 


^2 2 

*• 7 — cos k x cos k y + cos^ k x cos^ k y cos^ k z 

4 xs ys xs ys zs 


, D EH 4 , 2 , , DE^G 4 , 2 . 

+ — COS k x cos k z + — cos k x cos k y 

4 xs z s 4 xs x-’s 


4 2 2 

+ DEGH cos k x cos k y cos k z 
x s y's z s 


2 4 4 2 

+ DH G cos k x cos k z cos k y 
x s z s y^s 


2 4 4 

+ EG H cos k x cos k y 
x s y J s 


2 

cos 


k z 
z s 


(4.32) 
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The spatial average of equation (4.32) over the entire room volume is 
represented by an expression very similar to equation (4.14), That is, 


L L L 
x y z 


, 4 1 

% =• ■ V 


d> dx dy dz 
m s s s 


ooo 


Evaluation of this integral at the limits shown yields 


V> - its ° 2 * 2 - A » 2 » 2 + 


+ Tn g2h2 + 16 d2eh + Te de2g 


+ 32 DEG “ + m d “ 2g + 4s eg2h 


From Lyon's point source analysis. 


<^ m 4 > - (3/8) 3 . 


Therefore, in anticipation of comparing equation (4,34) to equation 

3 

(4.35), the quantity (3/8) is factored out of equation (4.34) to give 


<<j V 4> = + | d2 « 2 + § e 2 g 2 


+ g 2 h 2 + | d 2 eh + | de 2 g 


16 4 2 4 2 

+ -g DEGH + ^ DH G + ^ EG H 


(4.36) 
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Recalling once again the results of equations (4.18a) and (4.18b), 
it can be shown that 


lim 

d->o 


. 4 ,3 ' 

< 4 > > = («) 

m o 


9 


which equals the result shown in equation (4.35). Thus, as should now 

i | 

be expected, the quantity <<j> m > for the finite size model approaches 

4 

<ip > for a point source as the finite source dimension becomes very 
small. 

4 

Just as in the previous analysis of the mean, <<{> m > is a function 

of k , k , and k rather than the desired k . The substitution of 
x y z m 

equation (4.34) into equation (4.31), therefore, leads to an integral 

which is difficult to evaluate by direct means. Thus, the wave number 

4 

lattice is once again utilized to express <<p > as a function of k . 

m m 

This lattice average is given by 

tt/ 2 tr/2 

<4> k ^ sin 9 d 0 d y , (4.38) 

o o 



where the expressions developed in equations (4.22a), (4.22b), and 

(4.22c) have been substituted for k , k , and k . 

x’ y z 

The integration of this equation, which was once again attacked 
numerically using a digital computer, resulted in 



27 

256tt 


C(k d) 
m 


9 


(4.39) 
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where £ is understood to be a numerical value depending upon the magni- 
tude of k d. A brief description of the computer problem is given in 
m 

Appendix 9.1. 

Substitution of equation (4.39) into equation (4.31) yields the 
following expression for the variance of acoustic resistance for the 
finite source model: 



27 ,16pcgS? 
256-rr ' y 2 ; 


T k^n(k )£(k d)dk 

m m m 

[(k 2 -k 2 ) 2 + (^)Y 

o m V 


(4.40) 


Upon making the same assumption as in the previous analysis, i_.ji. 
the damping is small such that equation (4.40) peaks sharply about 
= k, and upon substituting 



and 


dR £ 2dk , 
m 


equation (4.40) transforms to 



2 2 

27 p c 
8 n „2 


n(k) E, (kd) 


(2gS/V) dR 

- [K 2 + C^)Y 


(4.41) 


Integration of equation (4.41), along with the utilization of equation 
(4,26) yields 



2 2 

27 p c 


16 i u 2 

kv n 


n(k) £(kd) 


(4.42) 
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Finally, normalization of equation (4.42) to the square of the mean, 
for a point source given by equation (3.30), leads to the expression 
for the normalized variance of power for the finite source model. The 
result is 


2 

„ 2 - (11 h 2 S( kd > 

2 N '“16 tt 

m„ 


(4.43) 


where M is the same modal overlap function as was defined in Chapter 3. 

It should be immediately evident that the quantity in parentheses 
in equation (4,43) is exactly the variance of power given in equation 
(3.4) for a point source. Thus, the variance of power for the finite 
size source can be put in the more illustrative form, 


where 


2 

°N (F.S.) 
2 

°N (F.S.) 
2 

°N (P.S.) 


2 

C N (P.S.) 


2g(kd) 

TT 


(4.44) 


variance of power for the finite source model, 


variance of power for a point source. 


If equation (4.44) is 
square root is taken, 
standard deviation of 
The result is 


2 

normalized with respect to cr . . and the 

JN * J e j 

the resulting equation expresses the ratio of the 
the finite source model to that of a point source. 


q N(F.S„) = f 2g(kd) 
a N(P.S.) V. 77 


(4.45) 



42 


which is plotted in Figure 4,5 using the values of £ determined from 
the computer data. The standard deviation ratio is always less than one 
and substantially less whenever kd is greater than or equal to three. 

In terms of acoustic wavelength, this value of kd corresponds to a 
finite source dimension greater than or equal to one-half an acoustic 
wavelength. Thus, the most significant result from Figure 4.5 is that 
the standard deviation for the finite source model is always less than 
the corresponding value for a point source and substantially less when- 
ever the ratio d/X is one-half or greater. 

In Figure 4,6, the theoretical standard deviation of power for 
several typical finite size sources is shown plotted against point 
source theory [equation (3,8)]. For simplicity in representation, the 
reverberation room is assumed to have a reverberation time of two seconds 
and a volume of 10,000 cubic feet. It can be seen that as the frequency 
or the source size increases, the normalized standard deviation of power 
decreases . 
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d/X 


4.5 Dependence of acoustic frequency or wavelength and source 

dimension on the ratio of the normalized standard deviation 
of power for a finite source to that for a point source 




NORMALIZED STANDARD DEVIATION OF POWER 
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FREQUENCY (Hz) 

4.6 Graph of the effect of source size on the uncertainty 
of pure tone sound power using a single source position 
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5 . EXPERIMENTAL INVESTIGATION 
5,1 Facilities and Test Setup 

The purpose of the experimental portion of this study was to gather 
data in support of the previous theoretical development for the finite 
size source model. Since a special reverberation facility was not 
available for use at the time of this study, the best hard room which 
could be found among the facilities of the Center for Acoustical Studies 
was selected for use. The dimensions of this room were 32.0 feet long 
by 26.3 feet wide by 9.6 feet high representing an internal volume of 
about 8,000 cubic feet. (Two column protrusions, representing a volume 
of 60 cubic feet, have been subtracted from the total volume.) These 
dimensions also correspond to a total room surface area for sound ab- 
sorption of about 2,900 square feet. Three of the room walls were 
painted cinder block, one wall was cinder block and glass, the floor 
was linoleum laid on concrete, and the ceiling was suspended acoustic 
tile with twelve large fluorescent light fixtures, A schematic diagram 
of this room is shown in Figure 5.1. 

One measure of the acoustical condition of the room is the average 
absorption coefficient of its walls. This parameter can be determined 
both theoretically using published absorption data and experimentally 
using Sabine’s diffuse room equation. Moreover, a good agreement 
between these two methods would be at least an indication of diffuse 
sound conditions. 

The theoretical determination of the average absorption coefficient 
is accomplished using the following well known equation (5.1): 
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5.1 Plan and elevation cutaway views of the reverberation facility 
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where a = the average absorption coefficient, 

= the absorption coefficient of the ith absorbing 
surface, 

= the surface area of the ith absorbing surface. 

The a^'s for use in equation (5,1) can be obtained from readily avail- 
able published data on the absorption coefficients of typical building 
materials. Table 5.1 shows the absorption coefficient and the associated 
area of each absorbing surface of the room as a function of the frequen- 
cies of interest in this study. The use of the data in this table along 
with equation (5.1) yields the following theoretical results: 


“125 

a 250 

a 500 


= 0.16, 
= 0.14, 
= 0.18, 


a 1000 = 0 . 22 , 


where the subscripts indicate the calculation frequency. These values 
indicate that this room would probably be classified as medium hard. 

The experimental average absorption coefficients can be calculated 
using the measured reverberation times of the room. Figure 5.2 shows 
the experimental setup which was used to make these reverberation 
measurements for one- third octave bands of random noise. The measured 
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Table 5.1 Absorption coefficients and surface area of reverberation 
room materials 


Absorption Coefficients Area 

125 250 500 1000 (ft 2 ) 


Acoustical tile 

0.44 

0.44 

0.61 

0.84 

660 

Cinder block 

0.10 

0.05 

0.06 

0.07 

1130 

Linoleum tile 

0.02 

0.03 

0.03 

0.03 

840 

Window glass 

0.09 

0.06 

0.04 

0.03 

90 

Light fixtures 

0.28 

0.22 

0.17 

0.09 

180 


Total surface area 


2900 




5.2 Test setup for reverberation time measurements 
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decay times can now be used in conjunction with Sabine's diffuse field 
relationship [equation (2.1)] to find the desired absorption values. 

Table 5.2 gives the results of the experimental study, again as a func- 
tion of the frequencies of interest. The theoretical and experimental 
values agree fairly well; suggesting the presence of a satisfactory 
diffuse field. 

As stated in Chapter 3, another measure of the room's applicability 
as a diffuse field is the Schroeder's large room frequency. Haying 
determined the average reverberation time of the room as 1.0 seconds, 
equation (2.2) yields 

f = 132 Hz . (5.2) 

By Schroeder's theory, a diffuse field exists above this frequency. 

The diffusivity of the room was checked in still another way using 
an eight inch speaker excited by random noise and placed in a corner of 
the room. Continuous sound level (A) data was taken of this random noise 
on two seven foot circular arcs at radial distances of ten and twenty 
feet from the source. The resulting graphs of sound levels, as shown in 
Figure 5.3, yield a difference in spatial average sound pressure of less 
than one decibel. This again indicates satisfactory diffuse sound 
conditions. 

The final room evaluation which needed to be made was a measurement 
of the ambient noise level. Figure 5.4 shows the one-third octave band 
frequency spectrum which was obtained. All band levels are relatively 
low except those in the range of 200 Hz and below. More will be said 
in a later paragraph concerning the precautions which were taken with 
respect to these relatively high levels in the lower frequency bands. 
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Table 5.2 Absorption coefficient values calculated using measured 
reverberation times and Sabine’s diffuse field equation 


Frequency 

(Hz) 

1/3 Octave 
Reverberation Time 
(sec) 

Absorption 

Coefficient 

125 

1.0 

0.13 

250 

1.0 

0.13 

500 

0.9 

0.15 

1000 

0.8 

0.17 

Average 

1.0 

0.14 




5.3 Continuous sound level (A) graphs at radial distances of 

10 and 20 ft from a random noise source in the reverberation 


room 


53 



5.4 Ambient noise spectrum in the reverberation room 
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The basic requirements for the experimental portion of this study 
included a system to generate the desired pure tones, a system to receive 
and analyze the resulting room pressure, and a device for recording the 
resulting data, Figure 5.5 is a schematic diagram of the complete equip- 
ment setup used to obtain the necessary data* The noise generation 
system consisted of a beat frequency or pure tone generator, a power 
amplifier, an ammeter, a frequency counter, and a sound source. The 
ammeter was used to keep the current output from the power amplifier 
constant for each test condition and, thus, regulate the source input 
power. 

The receiving and recording systems consisted of a microphone 
traversing apparatus, a spectrometer, and a graphic recorder. The 
traversing device, which was specially designed and constructed for use 
in this study, was used to facilitate the measurement of the spatial 
average sound pressure in the room, A complete description of this 
traverse is presented in Appendix 9.2, The second piece of equipment 
in this system, the spectrometer, served two purposes. First of all, 
it provided a means of measuring the actual sound pressure level detected 
by the microphone, and then it attenuated this microphone signal to a 
level acceptable for recording purposes. Secondly, through the use of 
its filtering system, the spectrometer eliminated the rather high com- 
ponents of low frequency ambient noise mentioned previously. The final 
instrument in the experimental setup was the graphic data recorder. 

This instrument provided a permanent record of the sound pressure fluc- 
tuations as detected by the traversing microphone. Thus, the data could 
be recorded and then analyzed at a later date. A photograph of the actual 
equipment setup is shown in Figure 5.6. 
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FREQUENCY 

COUNTER 


5.5 Experimental setup used to gather the statistical data 






56 



L 



57 


5.2 Test Procedure 

As was previously stated, the purpose of the experimental study 
was to obtain the data necessary to calculate the standard deviation 
of source power output from the mean output for low frequency, pure 
tone sound with source frequency, source size, and source location as 
the experimental parameters . As there is no instrument capable of 
measuring sound power directly, sound pressure must first be measured 
and then used to calculate the sound power. This relationship between 
sound pressure and sound power for any source located in a reverberant 
room is given by (4) 

PWL = SPL + 10 log 10 V - 10 log 1Q T - 29.9 , (5.3) 

-12 

where PWL = sound power level in dB, re 10 watts, 

SPL = spatial average sound pressure level in dB, 
re 0.0002 y bars, 

3 

V = room volume in ft , 

T = reverberation time in sec. 

Due to this proportionality of acoustic power to spatial average sound 
pressure, with all other quantities remaining constant for a given pure 
tone frequency, the variance associated with PWL is equal to that for 
SPL. Thus, to eliminate the chore of calculating the power associated 
with each spatial average datum value, only the statistics of sound 
pressure are actually calculated in the remainder of this study. 

It was decided to utilize four pure tone noise stimuli in the 
experimental study, namely 125 Hz, 250 Hz, 500 Hz, and 1000 Hz. It 
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should be noted that these selections correspond to the center frequen- 
cies of the third, fourth, and fifth standard octave bands, and, thus 
[from equation (3.6)], each successive frequency represents a fourfold 
increase in room modal density. The lowest excitation frequency, 125 Hz, 
was purposely chosen below the large room frequency of 134 Hz in order 
to test whether its corresponding statistics varied significantly from 
established theory for diffuse conditions. 

Two sound sources were selected for use, one an eight inch diameter 
speaker and the other a 30 inch diameter speaker, in anticipation that 
a change would be seen in the power statistics due to the difference in 
finite source size. Since both the point source theory and the finite 
source theory are based on a monopole type source, these speakers were 
baffled to best model this monopole generation, A photograph of the 
two sources utilized is shown in Figure 5.7. 

In selecting the source positions for use in this study, it was 
anticipated that only positions in which the speaker enclosure rested 
on the floor of the room would be utilized. In order for these source 
positions to be statistically independent, they had to be at least one- 
half a wavelength apart at the lowest frequency of interest. This, of 
course, amounted to 4.5 feet apart at 125 Hz. Also, due to the symmetry 
of the rectangular room, only source positions in one quadrant of floor 
area were considered. Using Beranek's (5) equation for the maximum 
number of source positions in a given space [modified slightly], 



(5.4) 
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5.7 Photograph of the eight inch speaker and the 30 inch 
speaker used in the experimental study 
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where N = maximum number of independent source positions, 

Q = area of floor quadrant, 

X = wavelength of lowest frequency, 

it can be determined that a maximum of eleven positions is possible. 
Obviously, the larger the number of source positions, the higher will 
be the confidence in the resulting data. Thus, the maximum number of 
source positions was utilized, and they were arbitrarily selected as 
shown by the room quadrant layout in Figure 5.8. 

In order to calculate the dependence of power statistics on source 
position for a given source, it is necessary to determine the spatial 
average sound pressure in the room for each source position. This 
spatial average is accomplished by measuring the sound pressure level 
at a number of independent locations in the sound field, i^. ei. , at points 
at least one-half a wavelength apart, and then obtaining the logarith- 
mic average of the resulting measurements. These microphone readings 
must be taken at a distance of at least one-half a wavelength from the 
source and at least one-fourth a wavelength from the room walls, since 
sound wave behavior in these areas deviates from the diffuse field 
condition. After careful consideration, it was decided that fifteen 
independent readings would provide a satisfactory approximation of the 
spatial average sound pressure level for each source position. The 
accuracy, of course, could have been improved with a larger number of 
microphone positions but the resulting increase in test effort by 
several orders of magnitude for just a small improvement in accuracy 
was not felt to be justifiable. 
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The required microphone readings were obtained as discrete indepen- 
dent points using the microphone traverse. Even with the use of this 
device, however, quite a bit of work was involved in that several 
different traverses were necessary for each source position. Recall 
from Chapter 2 that Lubman's expression for the maximum number of read- 
ings from a circular arc traverse is 


N = 1 + 2x/X , 


(5.5) 


where x is the traverse arc length. Using this relationship, the number 
of independent points for each frequency, given the traverse length to 
be 6.5 feet, is 


125 Hz 
250 Hz 
500 Hz 
1000 Hz 


2 points /traverse, 

3 points /traverse, 

7 points /traverse, 
12 points/traverse. 


The required number of traverses can now be found by dividing the desired 
number of independent readings [namely 15] by each of the above relation- 
ships. The result is 


125 Hz 

7.5 or 8 traverses , 

(5.6a) 

250 Hz 

- 5 traverses , 

(5.6b) 

500 Hz 

- 2.2 or 3 traverses , 

(5.6c) 

1000 Hz 

- 1.3 or 2 traverses 

(5.6d) 
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Two typical microphone traverse print-outs are shown in Figure 5.9. As 
is quickly noticed, fluctuations of 20 dB to 30 dB in the sound pressure 
level from point to point in the reverberant room are not uncommon. 

The actual data gathering process consisted of making fifteen sound 
pressure level readings for each source at eleven source positions and 
at four excitation frequencies. An attempt was made to regulate the 
speaker output in each test condition so that the average room sound 
pressure level was between 90 dB and 100 dB, a level at least 30 dB 
above the ambient noise. 

The final result of this test operation was the compilation of 
1155 data points, 660 points for the 30 inch speaker and 495 points for 
the eight inch speaker. 


5.3 Data Analysis and Results 

The compilation of experimental data is now analyzed using a 
standard two-way analysis of variance, namely an analysis both within 
and among source positions. This type of analysis can be facilitated 
by dividing the data into seven data point sets, with each set corre- 
sponding to a different combination of source size and excitation fre- 
quency. The result is 165 data points per set representing the different 
combinations of source position and microphone position. A typical data 
set, arranged in the form of a data point matrix, is shown in Table 5.3. 
Each element of the matrix y represents the squared sound pressure 
for the ith source location as measured at the j th microphone position. 


Due to insufficient speaker output, no data was taken for the 
eight inch speaker at 125 Hz. 
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1000 Hz 


5.9 Typical microphone traverse print-out data 
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Table 5.3 Data point matrix for a typical data set 
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Source 

Position 

1 

2 

Microphone 

3 

Position 

4 ... 

n 

1 

y ll 

y 12 

y 13 

y 14 

... 

y ln 

2 

y 21 

y 22 

y 23 

y 24 

... 

y 2n 

3 

y 31 

y 32 

y 33 

y 34 

... 

y 3n 

4 

y 41 

y 42 

y 43 

y 44 

... 

y 4n 

0 

. 

• 

• 

• 

• 

• 

• • 

• 

• 

• 

• 

• • • 

. 

• 

» 

• 

• 

• 

• 

• 

s 

y sl 

y s2 

y s3 

y s4 

• • • 

y sn 



The symbol n represents the total number of microphone locations and 
s denotes the total number of source positions. 
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Utilizing Table 5.3, two new quantities can be defined which will 
prove very helpful in mathematically describing the analysis of variance 
technique, These quantities are the spatial average sound pressure in 
the room for the source at the ith position, given by 



and the mean value of the spatial average pressures for all i source 
positions, given by 


s P 

G = l o“ • (5.8) 

i-1 


Consideration is first given to the average value of the normalized 
variance within source positions for a given data set. Based on a single 
microphone reading for each source position, this average variance can 
be expressed mathematically as 



1 

n-1 


n 

l 

j=l 


( y ir P i)‘ 


(5.9) 


This equation, which is applied to each data set, states, first of 
all, that the spatial average sound pressure in the room is calculated 
for the ith position of the source. The spatial average is then sub- 
tracted from each of the individual microphone readings y for that 
source position, and each result is then squared. The sum of these 
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squares over all n microphone positions, followed by division by 
2 

(n-1) , gives the normalized variance of a single microphone reading 

from the mean for the ith source position. If the above calculation is 

performed for all i source positions, and if the resulting variances 

2 

are averaged, then the average normalized variance a of sound pressure 
for a single microphone reading is obtained for that data set. A calcu- 
lation of this type is presented in Appendix 9.3 using an actual data 
set obtained from experimental investigation. 

The results of this analysis within source positions are given in 
Table 5.4 for each frequency and source size. It can be seen that both 
the variance and the standard deviation of the data are presented in 
both normalized form and in decibels. As to whether the normalized 
standard deviation or the standard deviation in decibels is more useful 
depends entirely upon the particular application. Thus, it is extremely 
helpful to be able to quickly find one form of the standard deviation 
given the other form. A graph showing the relationship between the two 
forms is given in Figure 5.10, 

It should be kept in mind that the results in Table 5.4 are based 
on a single microphone reading. As was mentioned in Chapter 2, if the 
spatial average sound pressure is based upon n microphone readings, the 
normalized variance is reduced by a factor of 1/n or 



2 . 
“a /n 


» 


(5.10) 



Table 5.4 Experimental uncertainty of pure tone mean square sound 
pressure using a single microphone position 
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Source 

Size 

(in) 

Frequency 

(Hz) 

Variance 

(dB) 

Standard 

Deviation 

(dB) 

Variance 

(normal) 

Standard 

Deviation 

(normal) 


250 

33.99 

5.83 

1.06 

1.03 

8 

500 

37.09 

6.09 

1.12 

1.06 


1000 

35.52 

5.96 

1.08 

1.04 


125 

28.20 

5.31 

0.92 

0.96 

30 

250 

36.48 

6.04 

1.10 

1.05 


500 

35.40 

5.95 

1.08 

1.04 


1000 


28.41 


5.33 


0.92 


0.96 
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5.10 Graph for converting normalized standard deviation 
to standard deviation in decibels 
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where 


2 

a = normalized variance due to a spatial average 
B 

based on n microphone readings, 

2 

a = normalized variance due to a spatial average 

A 

based on one microphone reading. 


The second part of the analysis of variance is performed to find 
the total variance or the variance among source positions. Based on 
n microphone positions per source position, this normalized variance 
of spatial average sound pressure [or sound power] from the mean pres- 
sure [or power] for all s source positions can be expressed mathemati- 
cally using Table 5.3, equation (5.7), and equation (5.8). The result 
is 


a 


2 

C 


1 

S-l 


s 

l 

i=l 


(VV 


(5.11) 


The mathematical operations involved in equation (5.11) are, first 

of all, the determination of the mean value G of the spatial average 

pressures in the room for all s source positions. This mean value is 

then subtracted in turn from each of the spatial average pressures P^. 

Each resultant is then squared, and a sum is taken over all s squares. 

2 

After division by G (S-l) , this sum represents the total normalized 
variance of source power from the mean based on the fact that the 
spatial average room pressure for each source position was determined 
using n microphone locations. 

As should be expected, therefore, the normalized variance given by 
equation (5,11) is actually the sum of two independent variances. One 
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variance is that due to incomplete spatial averaging of the resultant 
sound field for each source position. The value of this variance has 
already been determined and is given by equation (5.10). The second 
variance is that due to changes in source power caused by variations in 
source position. In other words, 





(5.12) 


where is the normalized variance due solely to changes in source 
position. The significance of this relationship may not be evident at 
this time, however in Chapter 6, where a comparison is made between 
theory and experiment, its utilization will be essential. 

A sample calculation of the total variance using an actual data set 
is given in Appendix 9.4. It should be emphasized again that the actual 
quantity determined in the Appendix is the variance of mean square pres- 
sure. However, due to the proportionality of mean square sound pressure 
to sound power, the quantity also represents the variance of sound power. 

The results of the statistical analysis of variance among source 
positions are given in Table 5.5. Again both the variance and standard 
deviation of each data set are presented in both normalized form and in 
decibels. 
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Table 5.5 Experimental uncertainty of pure tone sound power 
using a single source position 


Source 

Size 

(in) 

Frequency 

(Hz) 

Variance 

(dB) 

Standard 

Deviation 

(dB) 

Variance 

(normal) 

Standard 

Deviation 

(normal) 


250 

3 0 50 

1.87 

0.172 

0.415 

8 

500 

2.82 

1.68 

0.137 

0.370 


1000 

1.39 

1.18 

0.072 

0.268 


125 

5.38 

2.32 

0.260 

0.510 

30 

250 

0.92 

0.96 

0.048 

0.219 


500 

1.46 

1.21 

0.073 

0.271 


1000 

1.61 

1.27 

0.081 

0.284 
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6. COMPARISON OF THEORY AND EXPERIMENT 

Since the experimental results were presented in two parts - the 
variance within source positions and the variance among source positions, 
the comparison between theory and experiment will be presented in the 
same way. 

First consider the standard deviation of spatial average pressure 
within source positions for each data set. Recall from Chapter 2 that 
Doak in 1959 and Lubman in 1968 showed theoretically that this standard 
deviation for a pure tone source should equal 5.57 dB or, when normalized, 
should equal one. Comparing these theoretical findings to the experi- 
mental results of Table 5.4, the plots shown in Figure 6.1 for the eight 
inch speaker and in Figure 6.2 for the 30 inch speaker are obtained. For 
both source sizes, the agreement between theory and experiment is quite 
good. 

In order to test the hypothesis that the true finite source stan- 
dard deviation is equal to the theoretical value of one, a chi-squared, 

2 

X , distribution is assumed (10) . A test using such a distribution 
involves the definition of a new variable given by 


v 


(n-l)S 2 

2 


( 6 . 1 ) 


where v = experimental variable, 

(n-1) = degrees of freedom, 

S = experimentally determined normalized standard 
deviation, 


a 


= true standard deviation of the finite source. 
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0 0 



6.1 Uncertainty of pure tone sound pressure from an eight 
inch speaker using a single microphone position 


OQ 



6.2 Uncertainty of pure tone sound pressure from a 30 inch 
speaker using a single microphone position 
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For this case, the variable v possesses a x distribution with 164 

2 2 2 
degrees of freedom. If x^ and X£ represent the two values of x 

cutting off Z per cent tail areas, then the probability is (100 - 2Z) 

per cent that v will satisfy 


2 

X-l < v < 



( 6 . 2 ) 


Substituting equation (6.1) into equation (6.2), followed by algebraic 
manipulation, yields 


(n-l)S 2 2 (n-l)S 2 

2 < o < 2 • • 

x 2 x x 

2 

This expression is known as the Z per cent confidence interval for a . 
Taking the worst experimental case, namely S = 1.06 for the eight inch 
speaker at 500 Hz, the confidence interval becomes 


183.5 

2 


X 


2 


2 

< a 


183.5 



(6.4) 


2 

Arbitrarily choosing a confidence level of 95 per cent, the x tables 
indicate 


X ± 2 = 129.98 


(6 . 5a) 



200.86 


(6.5b) 


Substitution of equations (6.5a) and (6.5b) into equation (6.4), followed 
by the square root, yields 


0.96 < a < 1.19 


( 6 . 6 ) 



77 


for the 95 per cent confidence interval for the true finite source 

standard deviation [using the worst experimental case]. This interval 

includes the theoretical value of one, and, therefore the hypothesis 

is accepted with at least 95 per cent certainty. 

Secondly, the standard deviation of source power among source 

positions is considered. The normalized standard deviations of this 

type for a point source are determined theoretically using equation 

3 

(3.9). With a room volume of 8000 ft , this equation becomes 

2 T 1/2 

Vp.s.) * 3 ' 69 * 10 f • (6 ' 7> 

where the reverberation times are given in Table 5.2 for each frequency 

f. 

The theoretical standard deviations for both finite source sizes 

1/2 

are then found using equation (4,45). The values of (2£/tt) needed 
to evaluate this equation can be read from Figure 4,5, in which 

kd (30 inch speaker) = 1.40 x 10 ^f , 
kd (8 inch speaker) = 3.75 x 10 ^f 

The theoretical standard deviations determined for a point source 
and for the finite sources, however, are made assuming that the spatial 
average pressure in the room is known exactly for each possible source 
position. In other words, the theoretical values do not account for 
the incomplete spatial averaging inherent in the experimental results. 
Thus, before any comparison between theory and experiment can be 
attempted, the theoretical values must be adjusted to account for this 
extra experimental error. The required additive standard deviation due 
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to incomplete averaging can be found from equation (5.10) to be 

1/2 

= (1/15) = 0.258 , (6.8) 

c 

2 

or about 1„1 dB. Substitution of the theoretical values of (p g ) 

2 

and o . . and the value of c from equation (6.8) into equation 

N^r.boj d 

(5,12) then gives the adjusted theoretical values appropriate for com- 
parison with experimental results. 

The theoretical results for the two finite sources and for a point 
source are given in Table 6.1. This table includes both the adjusted 
and unadjusted values of standard deviation. [Recall that the experi- 
mental results for both finite sources was given in Table 5.5]. 

Figure 6.3 is a graph comparing the experimental and the adjusted 
theoretical results for the eight inch speaker to that predicted by 
point source theory. In addition to showing good agreement between 
finite source theory and experiment, both curves also fall near point 
source predictions. This should not be surprising due to the small 
value of kd for the eight inch speaker at low frequencies. 

The results for the 30 inch speaker are shown in Figure 6.4. 

Again, the agreement between the experimental results and the derived 
finite source theory is good. However, both finite source curves fall 
measurably below the standard deviations predicted by point source 
theory, especially at the two lower frequencies. The lack of a dramatic 
decrease at the two higher frequencies is due, at least in part, to the 
fact that the standard deviation attributable to incomplete spatial 
averaging [indicated by the dashed horizontal line in the figure] has 
become the dominant term in the total uncertainty. 
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Table 6.1 Theoretical uncertainty of pure tone sound power using a 

single source position (adjusted values are those including 
the error due to experimental incomplete spatial averaging) 


Source 

Frequency 

Standard 

Adjusted 

Standard 

Adjusted 

Size 


Deviation 

S tandard 

Deviation 

Standard 




Deviation 


Deviation 

(in) 

(Hz) 

(dB) 

(dB) 

(normal) 

(normal) 


125 

3.17 

3.50 

0.673 

0.720 


250 

1.53 

1.93 

0.337 

0.424 

point 

500 

0.67 

1.35 

0.152 

0.300 


1000 

0.37 

1.20 

0.073 

0.268 


125 

3.10 

3.43 

0.664 

0.712 


250 

1.45 

1.87 

0.320 

0.411 

8 

500 

0.52 

1.28 

0.124 

0.286 


1000 

0,10 

1.15 

0.030 

0.286 


125 

2.57 

2.85 

0.565 

0.621 


250 

0.67 

1.35 

0.155 

0.300 

30 

500 

0.12 

1.15 

0.032 

0.260 


1000 

0,05 

1.13 

0.002 

0.258 



NORMALIZED STANDARD DEVIATION OF POWER 


80 



FREQUENCY (Hz) 


6.3 Uncertainty of pure tone sound power from an eight inch 
speaker using a single source position 


STANDARD DEVIATION OF POWER (dB) 



NORMALIZED STANDARD DEVIATION OF POWER 
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6.4 Uncertainty of pure tone sound power from a 30 inch 
speaker using a single source position 


STANDARD DEVIATION OF POWER (dB) 
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Due to the limited number of source positions utilized in this 
study, it is important to determine the confidence which can be associ- 
ated with the conclusion that the standard deviation for the 30 inch 

2 

speaker is less than that for a point source. If a x distribution is 
again assumed with the variable v defined by equation (6.1), then the 
probability is (100 - Z) per cent that v will satisfy 

v < Xj 2 , (6.9) 

2 2 

where Xi is the value of x cutting off a Z per cent lower tail area. 

In other words, the upper confidence limit cutoff value for the true 
variance of the finite source data is 


2 . (n-l)S 2 
a < — 

X 1 


( 6 . 10 ) 


Taking the data point at 125 Hz where S = 0.621 as a typical value, then 
equation (6,10) with n = 11 source positions yields 


2 

a 



( 6 . 11 ) 


Assuming a confidence level of 70 per cent which corresponds to 
X x 2 = 7.50 , 

equation (6.11) yields 


a < 0.715 , 


( 6 . 12 ) 
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This upper limiting value is less than the standard deviation of 0.720 
predicted by point theory. Therefore, since the confidence level associ- 
ated with the 250 Hz point is much greater than 70 per cent, it may be 
stated with at least 70 per cent confidence that the standard deviation 
of power for the 30 inch speaker at 125 Hz and at 250 Hz is less than 
that predicted by point source theory. Unfortunately, the dominance of 
the uncertainty due to incomplete spatial averaging at 500 Hz and 1000 
Hz prevents the establishment of any meaningful conclusions from the 
experimental results at these frequencies. As previously stated in 
Section 5.2, the uncertainty due to incomplete spatial averaging could 
have been decreased with a larger number of microphone readings, but the 
huge increase in test effort necessary for a small increase in accuracy 
was not felt to be justified. 



84 


7. SUMMARY AND CONCLUSIONS 

Three significant conclusions can be drawn from the results of 
this study „ First of all, if two monopole type sources, one a point 
source and the other the finite source model, have equal volume veloci- 
ties, their diffuse field acoustic power output is also the same, 
Furthermore, the diffuse field power outputs for both source sizes are 
equal to the corresponding powers as measured in the free field. 

Secondly, for a finite source, the standard deviation from the 
mean of pure tone sound pressure at a point is equal to that derived 
for any pure tone source located at some given source position. Further- 
more, this standard deviation is independent of source size or source 
frequency and has a decibel value of 5.57 dB or a normalized value of 
one. 

Finally, the normalized standard deviation of low frequency, pure 
tone sound power output for a monopole type finite size source is lower 
than that predicted for a point source. This decrease is considerable 
when the ratio of source dimension to acoustic wavelength is one-half 
or greater. The significance of this conclusion is that in some cases 
the reliable measurement of sound power output (within about - 2 dB) 
of a low frequency, pure tone finite size source can be made without 
utilizing the large number of source positions required by point source 
theory. In other words, the use of a reverberation room to measure the 
power output of a finite size source, even for source spectrums contain- 
ing a large, low frequency pure tone or narrow-band component, may prove 
to be quite attractive. This should especially prove to be true for 
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rooms equipped with rotating diffusers which decrease the sound pressure 
fluctuations caused by the room modes , 

In addition to the above conclusions, another contribution of this 
study is the development of a microphone traversing apparatus to facili- 
tate the determination of spatial average pressure. This traversing 
device was used extensively in this study and could prove quite useful 
in future reverberant field research studies at the Center for Acoustical 
Studies, To the author's knowledge, such a traversing mechanism is not 
available in the commercial market. 

While this study has provided understanding of one of the important 
missing links in the use of reverberant rooms for sound power measure- 
ments, additional extensive research is needed before a satisfactory 
"total" understanding is achieved. For example, the work done in this 
study and the majority of work done concerning point sources has been 
established assuming the source to have monopole type radiation charac- 
teristics. The statistics of sources having dipole, quadrapole, and 
higher order characteristics are essentially unknown. However, two 
research projects concentrating on dipole behavior are presently under- 
way. Although the results are not yet available, Yousri and Fahy at 
the University of Southampton are studying the statistics of point dipole 
sources in reverberant rooms. In another effort, Lubman is studying the 
behavior of the ILG or "random dipole" source. Although once again the 
results are not yet available, the claims that preliminary indications 
show a higher standard deviation than for a point monopole source, at 
least at low frequencies. Both of these research efforts should provide 
still another stepping stone toward the establishment of a high level of 
confidence in reverberant room measurements. 
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In addition, as far as finite size sources are concerned, this 
study investigated only one type of source orientation - that parallel 
to the y-z plane. The statistics of other orientations, such as a skew 
orientation, have not been investigated. Finally, Eldred^ has suggested 
that the modal density of reverberant rooms excited by large sources 
might be decreased by several orders of magnitude at high frequencies 
due to poor coupling between the source and the sound field. Although 
this theory has not been fully substantiated, it could foreseeably add 
another factor of confusion into the many problems already known. 

However, even with the present existing problems, it would appear 
that in the near future the reverberant room may emerge as the most 
extensively utilized device for any type of acoustical measurement not 
requiring directivity information. 


The concept of decreased modal density was brought out in "Effective 
Modal Density in a Reverberant Sound Field for Finite Sized Sources" by 
B. Ho Sharp and K. M, Eldrea, presented to the Acoustical Society of 
America, April 18, 1972. 
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9. APPENDICES 

9.1 Numerical Integration Procedure 

The objective of the numerical integration procedure is to obtain 
good approximations for equation (4.24) and for equation (4.28) as a 
function of the quantity k^d. Due to the similarity of these two 
equations, the same numerical procedure is utilized to integrate both 
expressions. Therefore, the actual computer programs are essentially 
alike for both approximations, except for the obvious difference in 
the two integrands. 

The procedure developed consists of approximating each integral 
as a finite double sum over 0 and y with the limits for both sums being 
from zero to -rr/2. The accuracy of each approximation is controlled by 
the number of point values included within each summing interval. 

For simplicity in describing the computer programs, let k represent 

2 4 

both «{> > sin 0 in equation (4.24) and <<j> > sin 0 in equation (4.28), 

where it is understood that k is a function of 0, y, and k d for both 

m 

cases. A conventional digital computer flow diagram can then be used 
to describe the basic steps in the numerical procedure. This flow 
diagram is shown in Figure 9.1. 

9.2 Microphone Traverse 

Even before this study was undertaken, it was determined that a 
microphone traverse mechanism was needed to facilitate the determination 
of spatial average pressure in the reverberation room. It was intended 
that this device would move a microphone either through a straight line 
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or through a circular arc, and thus provide several independent pressure 
readings at one time. 

In the preliminary stages of design, several criteria were estab- 
lished to ensure that the device would meet expectations. First of all, 
the need for a means of adjusting the traverse length was felt to be 
desirable. The use of this device in different sized rooms and for 
different frequency ranges might warrant such an adjustment. Secondly, 
the traverse had to make a minimum of noise during operation. Any opera- 
tional noise could interfere with the acoustical signal being detected. 
This criterion placed quite a restriction on the design, since it prac- 
tically prohibited the use of any motorized device. Next, the device 
was to be controllable from outside the reverberation room. It was 
intended that this stipulation would prevent any operator disturbance 
of the sound field by enabling the operator to start and stop the device 
at will without having to enter the reverberation facility. It was also 
felt that the microphone motion itself should be as constant as possible 
so that repeated tests under the same conditions would produce the same 
data print-outs. Next, economic factors had to be considered, since 
funds were not available for a large expenditure. The design had to be 
geared toward construction using existing materials or ones which were 
inexpensive and readily available. Finally, the traverse had to be 
durable enough to withstand continual movement from one room location 
to another and yet be light enough to facilitate relocations without 
much effort. 

Before attempting to completely design such a device, a search was 
made to see if a device meeting these specifications already existed at 
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any of the well known research centers in the country. One such center 
claimed to have such a device, and, therefore, a complete set of working 
drawings was purchased from them. These plans were quite helpful in 
providing a good design basis, but as far as a practical and usable 
final design was concerned, the plans were rather incomplete and the 
finished produce highly unacceptable. Many design alterations were 
made before an acceptable traversing apparatus was obtained. 

While the original plans called for a hydraulic device actuated by 
pressurized air, the device actually constructed operates entirely on 
pressurized air. Boom vertical motion is produced by the continuous 
compression and decompression of a standard type air cylinder. Control 
of the motion is accomplished using mercury switches on the traverse 
boom and an electrical relay in the control box which actuates the air 
inlet and exhaust valves. Since both of the valves and the relay are 
located outside the test room, the only noise produced by the traverse 
mechanism itself is the frictional noise of the air cylinder piston 
and of the boom bearing. 

The length of the boom itself provides the necessary control over 
the microphone sweep length. Furthermore, this boom as well as the 
majority of the device, is constructed of aluminum, thereby limiting 
the total weight to about twenty pounds. Most of this weight is 
attributable to two inertial masses located on the boom which serve to 
provide smooth motion at a moderate angular velocity. 

Figure 9.2 shows a simplified drawing of the basic components of 
the microphone traverse device. This is followed in Figure 9.3 by a 
drawing of the unit designed to control the traverse motion. 
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9.2 Simplified drawing of the microphone traverse apparatus 


vo 



ON-OFF 

SWITCH 


-AIR EXHAUST VELOCITY CONTROL 

^^AIR HOSE TO TRAVERSE 



TRAVERSE 
MOTION 
LIGHTS 


cr 


i 


I 






UP 

DOWN 

RELAY 

CONTROL 

BOX 


AIR INLET 
VELOCITY CONTROL 


AIR 

EXHAUST 
VALVEy 


AIR 

INLET 

VALVE 




H==L 



PRESSURE REGULATOR (30 PSIG) 
TO MERCURY SWITCHES 



9.3 Simplified drawing of the control unit 





The specifications of the traverse apparatus as it was finally 
constructed can now be determined utilizing Figure 9,4 which shows the 
traverse boom motion, greatly simplified. In list form, the specifica- 
tions are: 


Traverse arc angle = 53°, 

Boom length = 7 ft, 

Traverse arc length = 6.5 ft, 

Boom velocity = 20 traverses in 35 sec 
= 1,75 sec/traverse, 

Delay time for boom direction change z 0.10 sec. 

Time for one traverse = 1.75 - 0,10 = 1.65 sec. 

Boom velocity =6.5 ft/1.65 sec = 4 ft/sec. 

In terms of a graphic recorder print-out with a paper speed of 30 mm/sec, 
this boom velocity corresponds to a traverse sweep length on paper of 
approximately 2.06 inches. 

It should be understood that the above specifications are valid 
only when the traverse device is equipped with the seven foot boom 
utilized throughout this study. Should this boom length be increased 
or decreased, the specifications would change accordingly, 

A photograph of the actual traverse device and its control unit are 
shown in Figure 9.5. As can be seen, both components of the total system 
are constructed of parts which are relatively inexpensive and easily 


obtained. 
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9.4 Simplified drawing of microphone boom motion 
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9.5 Photograph of microphone traverse mechanism and its control unit 
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9.3 Sample Calculation: Standard Deviation 

Within Source Positions 

The procedure for calculating the standard deviation within source 
positions for a given data set was explained in detail in Chapter 5. 
Table 9.1 shows actual data taken for the 30 inch speaker at 250 Hz and 
at source position seven. Included in the table are the necessary 
mathematical calculations involved in the computation of the required 
standard deviation. The following symbols are used in the table. 


SPL. 

i 


= sound pressure level at ith microphone location, 



( p/p r )2 

SPL 


= ratio of the actual acoustic pressure at microphone 
position i to the reference pressure of 0.0002 
dy/cm, 

= spatial average squared sound pressure, 

= spatial average sound pressure level. 



= variance and standard deviation in dB, respectively, 

= variance and standard deviation normalized to the 
mean square and the mean, respectively. 


Having completed the computation shown in Table 9. 1 for source 
position seven, the same procedure is utilized to analyze the remaining 
ten source positions. All eleven of the resulting uncertainties are then 
averaged to obtain the mean standard deviation within source positions 
for this data set. Table 9.2 shows the actual data, the calculated values, 
and the resulting mean standard deviation for the 30 inch speaker at 250 Hz. 
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Table 9.1 Sample calculation of the uncertainty of mean square 

pressure for a single microphone position using source 
position seven for the 30 inch speaker at 250 Hz 


Microphone 

Position 

SPL ± 

(SPL-SPL) 1 

xlO 8 

t(p/p R ) "(P/P 
*10 16 

1 

79 

121.0 

0.80 

82.68 

2 

85 

25.0 

3.16 

45.32 

3 

89 

1.0 

7.95 

3.77 

4 

93 

9.0 

20.00 

102.16 

5 

89 

1.0 

7.95 

3.77 

6 

89 

1.0 

7.95 

3.77 

7 

81 

81.0 

1.26 

74.52 

8 

85 

25.0 

3.16 

45.32 

9 

94 

16.0 

25.18 

233.70 

10 

95 

25.0 

31.60 

471.21 

11 

85 

25.0 

3.16 

45.32 

12 

94 

16.0 

25.18 

233.70 

13 

87 

9.0 

5.02 

23.74 

14 

80 

100.0 

1.00 

79.08 

15 

87 

9.0 

5.02 

23.74 


100 


Table 9.1 (Continued) 


l ( P /p_) 2 xl0 8 /15 = 9.89xl0 8 = (p/ Pp ) xlO 8 
• 1 ^ 1 *'■ 


f " Q . _ . 

10 log (p/p B ) xlO = 89.97 dB = SPL 
K 

°dB = £ (SPL-SPL) 2 /14 = 33.14 dB 2 

a . £=1 


adB 


to 2 y 


1/2 


= 5.75 dB 


<4 " l [(p/P R ) 2 “(p/P R ) 2 ]i xlO 16 /[14x(p/p r )" 4 xlO 16 ] = 1.07 
- 1/2 

°AN " “ M 3 
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Table 9.2 Sample results for the average standard deviation of mean 
square pressure for a single microphone position using 
the 30 inch speaker at 250 Hz 


Source 

Position 

Standard 

Deviation 

(dB) 

Standard 

Deviation 

(normal) 

1 

7.81 

1.22 

2 

5.40 

0.98 

3 

5.98 

1.04 

4 

8.24 

1.26 

5 

4.24 

0.83 

6 

7.00 

1.14 

7 

5.85 

1.03 

8 

6,06 

1.05 

9 

5.02 

0.93 

10 

6.82 

1.12 

11 

4.08 

0.82 

Average 

6.04 

1.04 


9.4 Sample Calculation: Standard Deviation 

Among Source Positions 
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The procedure for calculating the standard deviation among source 
positions was given in detail in Chapter 5. Table 9.3 shows the spatial 
average sound pressure levels and the spatial average squared sound 
pressures which were calculated from data obtained for the 30 inch 
speaker at 250 Hz. Also included are the mathematical calculations 
which are necessary to yield the required standard deviations. The 
following symbols are used in the table. 


(p/Pr ) 2 

3 


= spatial average squared pressure for the jth source 
position, 


SPL. 

1 


= spatial average sound pressure level for source 
position j , 


SPL 



= mean value of the spatial average squared pressures, 

= mean value of the spatial average sound pressure 
levels , 

= variance and standard deviation in dB, respectively, 
= variance and standard deviation normalized to the 


square of the mean and the mean, respectively. 
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Table 9o3 Sample calculation of the uncertainty of sound power for a 
single source position using data taken for the 30 inch 
speaker at 250 Hz 


Source 

Position 


(SPL-SPL) 2 

_ _ _ O 

O — n/1 O 

SFL j 

( p/p R )j 

[( p /p R )‘'-( p /p R )‘']“ 

1 

86,75 

4,750 

4,74 

9,550 

2 

88,78 

0,020 

7.55 

0,078 

3 

88,05 

0.770 

6,41 

2,016 

4 

88,06 

0,176 

6,42 

1,988 

5 

89,35 

0.176 

8,63 

0.640 

6 

89.42 

0,240 

8,75 

0.846 

7 

89,97 

1,080 

9.89 

4.243 

8 

88,89 

0,002 

7,75 

0,006 

9 

89.80 

0.757 

9.57 

3,028 

10 

88,20 

0.532 

6,61 

1,488 

11 

89.92 

0,980 

9.83 

4,000 


l (p/Pp) j x 10 8 /11 = 7 o 83x10 8 = (p/p ) 2 xl0 8 
j-1 J R 

10 log(p/p R ) 2 xl0 8 - 88 o 9 3 SPL 

2 l 1 __ 2 o 

o dB =1 (SPL-SPL) ^ /10 = 0.950 dB Z 
j=l J 
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Table 9 0 3 (Continued) 


a 


dB 


< 4 > 


1/2 


0 „974 dB 


°CN = ^ [ (P /P R ) 2 ~ (p /P R )^ ] 2 xl0 16 / [ lOx (p/P R ) ^xlO 16 ] - Oo 045 

j=l 

2 1/2 


0 o 213 



